Imprint of the merger and ring-down on the gravitational wave background from 

black hole binaries coalescence 
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We compute the gravitational wave background (GWB) generated by a cosmological population 
of (BH-BH) binaries using hybrid waveforms recently produced by numerical simulations of (BH- 
BH) coalescence, which include the inspiral, merger and ring-down contributions. A large sample of 
binary systems is simulated using the population synthesis code SeBa, and we extract fundamental 
statistical information on (BH-BH) physical parameters (primary and secondary BH masses, orbital 
separations and eccentricities, formation and merger timescales). We then derive the binary birth 
and merger rates using the theoretical cosmic star formation history obtained from a numerical study 
which reproduces the available observational data at redshifts z < 8. We evaluate the contributions 
of the inspiral, merger and ring-down signals to the GWB, and discuss how these depend on the 
parameters which critically affect the number of coalescing (BH-BH) systems. 

We find that Advanced LIGO/ Virgo have a chance to detect the GWB signal from the inspiral 
phase with a (S/N) — 10 only for the most optimistic model, which predicts the highest local 
merger rate of 0.85 Mpc"'' Myr~^. Third generation detectors, such as ET, could reveal the GWB 
from the inspiral phase predicted by any of the considered models. In addition, ET could sample 
the merger phase of the evolution at least for models which predict local merger rates between 
[0.053-0.85] Mpc"^ Myr'\ which are more than a factor 2 lower the the upper limit inferred from 
the analysis of the LIGO S5 run [J] . 

The frequency dependence and amplitude of the GWB generated during the coalescence is very 
sensitive to the adopted core mass threshold for BH formation. This opens up the possibility to 
better understand the final stages of the evolution of massive stellar binaries using observational 
constraints on the associated gravitational wave emission. 

PACS numbers: 04.30.db, 04.25.dg 
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I. INTRODUCTION 



Double black hole binaries are among the most promis- 
ing sources of gravitational radiation for the ground- 
based detectors LIGO/ Virgo - in their present and ad- 
vanced configurations which plan to increase the sen- 
sitivity by a factor 10 (http: / /www.ligo.caltech.edu/[ 
http://www.ego-gw.it/) - and for the future interferom - 
ctcrs, the space detector LISA (http://sci.esa.int/lisa) 
and the Einstein Telescope (ET, www.ct-gw.eu). In ad- 
dition to the emission from single resolved systems, mas- 
sive compact binaries generate stochastic backgrounds of 
gravitational waves (GWB), as extensively discussed in 
the hterature ■ If detected, these backgrounds would 
provide information on the cosmic star formation history, 
on the evolution of compact stars in binary systems, and 
on cosmological parameters [l3, El ; even if these signals 
were not detected, more stringent upper limits would rule 
out the most optimistic theoretical models. Astrophys- 
ical backgrounds act as a confusion-limited foreground 
noise for signals generated in the early Universe (pri- 
mordial GWBs); therefore, the spectral properties which 
characterize the background generated by various fami- 
lies of compact binaries need to be accurately modeled, 
and specific techniques need to be envisaged in order 
to disentangle their contribution from the instrumental 
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The aim of the present work is to provide updated es- 
timates of the GWB generated by a cosmological popula- 
tion of (BH-BH) binaries, including the contribution due 
to the merging of the two bodies, and to the ring-down 
of the final black hole. The full signal is modeled using 
the waveforms derived in [T5| : they combine the inspi- 
ral part of the signal (obtained with the standard Post- 
Newtonian description) and the ring-down oscillations of 
the final black hole, with the signal emitted during the 
merging phase, whose description has been made possi- 
ble by recent progresses in numerical relativity. We will 
show that since a significant amount of energy is radi- 
ated during the merger and ring-down, these phases give 
a significant contribution to the produced background. 

A recent estimate of GWB from coalescing black hole 
binaries using the same waveforms has been presented 
in 0. This study adopts average quantities for the sin- 
gle source emission and a fixed distribution for the delay 
between the birth of a binary and its merging (merger 
timescale); furthermore, they normalize the (BH-BH) 
merger rate to a "local value" inferred from models of 
Galactic binary populations. 

Here we follow a different approach. We 
use the population synthesis code SeBa [Tej 
( http://www.sns.ias. edu/~ starlab/| to simulate the 
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properties of a large sample of double black hole binaries. 
This enables us to investigate the statistical distributions 
of masses and merger timescales in a self-consistent way. 
Following [13, [ill , we adopt a theoretical model for the 
cosmic star formation history at redshifts z < 20 taken 
from the numerical simulation of [l9l |. which reproduces 
the observational data available at z < 8 [20]. With 
these inputs, we compute the redshift evolution of the 
(BH-BH) birth and merger rates and the cumulative 
gravitational wave signal produced throughout the Uni- 
verse. We evidenciate the contributions of the inspiral, 
merger and ringdown phases of the coalescence process 
to the total background, and discuss their detectability. 

Finally, we run a set of independent SeBa simulations 
to single out the parameters which critically affect the 
number of coalescing (BH-BH) systems, the GWB and its 
detectability; we find these to be the common envelope 
parameter, the core mass threshold for BH formation and 
the kick velocity distribution. 

As first noted in [2i|, close black hole binaries form 
efficiently through dynamical interactions in dense glob- 
ular clusters or in nuclear star clusters [2^ - [2^ . However, 
in the present analysis we do not consider this additional 
formation channel, since it requires a different modeling 
which is not included in the current version of SeBa. 

The plan of the paper is the following. In section [IT] 
we briefly describe the SeBa code and the underlying 
physical assumptions which we use to simulate the sam- 
ple of (BH-BH) binaries. A statistical analysis of the 
resulting population is presented in section IHII In sec- 
tion IIVI we outline the main features of the waveforms 
we adopt to describe the gravitational emission of (BH- 
BH) binaries and in section fVl we compute, starting from 
the cosmic star formation history, the evolution of the 
birth and merger rates of (BH-BH) binaries. In section 
IVIl we present the resulting density parameter of the 
GWB, JIqw, we discuss its detectability by second and 
third generations interferometric detectors, and its de- 
pendence on key physical parameters. Finally, in section 
IVIII we draw our conclusions. 



II. SEBA POPULATION SYNTHESIS CODE 

To compute the statistical properties of a black hole 
binary population we use the latest release of the popula- 
tion synthesis code SeB a Il6l. w hich is based on previous 
versions described in [3, 125l - l27l | . 

In SeBa. the evolution of binary systems is modeled by 
taking into account the relevant physics, which include 
stellar composition, stellar wind, mass transfer and accre- 
tion, gravitational radiation, magnetic braking, common 
envelope phase and supernovae (for more details see the 
original papers). The present version of the code uses up- 
dated stellar and binary physics, including results from 
supernova simulations [28|. In particular, new features 
about the evolution of massive stars, Wolf-Rayet stars, 
stellar wind mass loss rate, common envelope phase, fall- 



back prescription and supernova kicks are discussed in 
detail in [ig. The predictions of the code have been 
tested against observational data on double neutron star 
binaries [26|, BH-candidate short-period binaries [l6j . 
Type II and Ib/c SNe ll,!!^, and more recently, SN Type 
la [23. Here we sketch a summary of the assumptions 
that are relevant for the evolution of (BH-BH) binaries. 

For the common envelope evolution (CE), we use the 
standard prescription described in [3, 113, HH, with the 
efficiency and the structure parameters acE and A com- 
bined into a single quantity, acE x ^- The treatment 
of the common envelope evolution is still under debate; 
at present, a strict criterion to define the binding energy 
of the stellar envelope is still lacking and it is unclear 
whether other sources of energy, beyond gravitational en- 
ergy, contribute to unbind the common envelope ( [s^ - iSTj 
and references therein). As a result, acE x A is a free 
parameter of population synthesis models. For our refer- 
ence model A (see Section lHl]) . we fix acE x A = 2, chosen 
so as to reproduce the available observational constraints 
jssj , such as the (NS-NS) merger rate inferred from dou- 
ble pulsar observations |39l . |40[ . 

As suggested in [iTMisj . black holes receive a small 
asymmetric kick at birth: in our reference simulation 
(model A) this kick is taken from a Paczyiisky velocity 
distribution [i^, , isotropic in space and scaled down 
with the ratio of black hole to neutron star masses (see 
Table U). 

Black hole formation is treated in the code assuming 
that a constant fraction of the supernova explosion en- 
ergy is used to unbind the stellar envelope ; we choose 
/ = 0.4. However, while in the explosion energy is 
assumed to be a function of the pre-supernova mass, we 
keep it fixed at 10^° erg. This value is within the expected 
range, but favours the formation of rather massive black- 
holes (up to 15 AIq, see also [l6jV 

Regarding the evolution of massive (15 Mq < Mi < 
85AIq) stars, SeBa assumes that stellar wind mass-loss 
rate increases in time; in their total lifetime, stars lose 
an amount of matter which is a function of their initial 
mass, O.OIM^; the hydrogen envelope is lost when the 
stars are still on the main sequence. For higher masses 
{Mi > 85Mq) the mass loss rate is very uncertain: we 
assume that, during the main sequence lifetime, these 
stars lose 43 Mq. Due to the paucity of these massive 
stars, this crude assumption has negligible consequences 
on the estimate of the GWB. 

When the star loses its hydrogen envelope on the main 
sequence, stellar mass loss prescription for Wolf-Rayet 
stars is adopted. More details on the treatment of mas- 
sive stars in the code are described in [l^ (see in par- 
ticular section 2.3 and Fig.l). For Wolf-Rayet stars the 
mass loss rate by |4^ is adopted, which is based on the 
compilation of mass-loss rates inferred from observations 
of Wolf-Rayet stars [47 1. 

Wc initialise iV = 10'^ "zero-age" binaries (ZAMS). 
The zero-age parameters of the simulated population are 
given in Table |T] for our reference model A. These are 
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TABLE I: Zero-age parameters for the reference Model A (see text). 



Model A 


Parameter 


Symbol 


Value 


Note 


Mass of primary star 


.^^prim 


[ 8-100] A/© 


Salpeter IMF (-2.35) 


Mass of the secondary star 


Msec 




the distribution matches the q-distribution 


Mass ratio 


q 


[0-1] 


flat distribution 


Initial semi-major axis 


sma 


0.1 - 


flat distribution in log sma 


Eccentricity 


e 


[0-1] 


thermal equilibrium distribution 


CE parameter 




2 


structure parameter 


Kick distribution for NS 


u — vja 


a = 300 km s"^ 


Paczyriski distribution for « 


Kick distribution for BH 


"BH 




same as NS but scaled down: ubh = ''(Mns/A^bh) 


Core mass threshold 


Wlthro.BH 


IOM0 


fl6] 



randomly selected from a set of independent distribution 
functions. In particular, the initial primary mass Afpi-im 
is taken from a Salpeter Initial Mass Function (IMF), 
$(M) oc M-(i+^) with X = 1.35, between [8-100]!/©; 
the initial secondary mass, Mgoc, is selected from a flat 
distribution for the mass ratio q = Afsec/Afprim- The 
semi-major axis (sma) distribution is flat in log(sma) [48| 
ranging from 0.1 Rq (Roche lobe contact) up to lO^i?©. 
We assume a thermal eccentricity distribution $(e) = 2e 
in the range [0-1] [i^. Kicks follow a Paczynski velocity 
distribution with a dispersion of cr = 300 km s^^ [iij . 

We have explored a wide parameter space to investi- 
gate how different choices affect the resulting GWB. This 
has enabled us to identify the physical parameters which 
have the largest impact on the GWB, and the correspond- 
ing models will be presented in section fVIB I 

Finally, we note that the version of SeBa that we use 
assumes an initial solar metallicity for all stars. 



III. STATISTICAL PROPERTIES OF BLACK 
HOLE BINARIES POPULATION 




■ coalescing systems 
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FIG. 1: Mass of the secondary stellar progenitors as a func- 
tion of the corresponding primary mass. Empty grey circles 
indicate stellar binaries which lead to (BH-BH) systems; black 
squares show the subsample of these progenitors which form 
(BH-BH) binaries with merger times smaller than the Hubble 
time (see text). 



In this section, we present the statistical properties 
of (BH-BH) binaries and their progenitors for our ref- 
erence model A (see Table U) . A great uncertainty in 
the modeling of black hole binaries is due to our poor 
knowledge of the BH mass distribution and its relation 
to the initial distribution of progenitor masses. The up- 
per mass limit of isolated stars for black hole formation 
depends predominantly on wind mass loss; recent stud- 
ies constrain the mass range of black hole progenitors to 
M ~ [2O-6O]M0 [2l[3i,[5fl[5l|. Clearly, a massive pro- 
genitor in a binary system might follow a different evolu- 
tionary path which would affect the mass of the nascent 
black hole H El III HI . 

Figure [1] shows the mass range of primary and sec- 
ondary progenitors of (BH-BH) systems (empty grey 
circles). Black squares show the sub-sample of pro- 
genitors which generate (BH-BH) binaries with merger 
times smaller than the Hubble time (^h ~ 13.4 Gyr 
for our adopted cosmological model [Z^). Primary (sec- 
ondary) progenitors have masses in the range 30Mq < 
Mprim < IOOA/q (2OA/0 < Msec < IOOMq). Merging 



systems come from a narrower dynamical range, with 
90% (96 %) having primary (secondary) progenitor mass 
30Mq < Mprim < 6OM0 (20Mq < Msec < 5OM0). 

Observations of BH candidates in binary systems sug- 
gest a range of masses in the [4 — 17]M0 interval 
([H, Ulillll and references therein), in broad agreement 
with theoretical simulations, which indicate a range of 
black hole masses up to [10 - 15]M0 [H, [H]. In Fig. H 
we plot the primary and secondary black hole masses, 
Afi^BH and M2.BH (black squares) for model A. Primary 
and secondary black holes have comparable masses, rang- 
ing between ~ 6M0 up to ~ 2OM0; these limits are con- 
sistent with the observational and theoretical estimates 
quoted above. 

Merging pairs are concentrated in a narrow range of 
masses, with the primary BH mass in the range [10 — 
15]M0 and the secondary in the range [10 — 11]M0. This 
corresponds to a BH mass ratio, geH, in the range [1 - 
1.1] (see also Fig.|4]), where ^bh is defined as the ratio of 
the heaviest to the lightest BH mass. 

In the upper panel of Fig. [31 wc show the time interval 



4 




FIG. 2: Same as Fig[l] but for black hole primary and sec- 
ondary masses. The various regions in the parameter space 
which are over- or under-populated originate from the deci- 
sion making process in SeBa, and can be related to the chain 
of events in the population synthesis. 



from the formation of the ZAMS binary system to the 
formation of the compact black hole binary, Tg; we plot 
this quantity as a function of the semi-major axis, sma, 
showing only those systems with sma < 1600 Rq, which 
represent ~ 60% of the total. The formation time of 
compact binaries is very small, ranging between ~ 3.5 to 
~ 6 Myr. Black squares indicate merging black hole bina- 
ries, which arc those that, at the time of their formation, 
have semi- major axis smaller than 20 Rq , comparable to 
what found in (see their Fig. 2). 

It is important to note that only 2% of (BH-BH) pairs 
is able to reach the final coalescence. Figure [3] shows 
that the majority of (BH-BH) systems is characterized 
by large orbital separations (and periods). In fact, pro- 
genitors with masses > 40AfQ experience large mass loss 
rates which remove mass from the binary increasing the 
orbital separation. 

Once compact degenerate binaries are formed, the 
emission of gravitational radiation is the only physical 
process driving the change in orbital parameters. In the 
lower panel of Fig. [3l we show the merger time as a func- 
tion of the chirp mass, defined as 

where ^ = Mi^bhM2,bh/M and M = Mi^bh + Af2,BH 
arc the reduced and total mass, respectively. It is clear 
from the figure that all merging pairs have relatively long 
merger times, > 1 Gyr and 70% have Tm > 6 Gyr (hor- 
izontal dashed line); the chirp masses lie in the interval 
[8.5- 1O.5]M0 and 75% of the systems have M < 9.6Mq 
(vertical dashed line). Although the plot does not evi- 
denciate any clear correlation between Tm and M, half 
of the systems lie in the upper left region and only 5% in 
the bottom right. 

Finally, in Fig. 2] we show the distribution of orbital 
eccentricity as a function of the black hole mass ratio for 
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FIG. 3: Upper panel: formation timescales for compact black 
hole binaries as a function of the corresponding semi-major 
axis. For clarity, only systems with sma < 16007?© are shown 
(empty grey circles). Black squares are (BH-BH) pairs with 
merger times smaller than the Hubble time (see text). Lower 
panel: merger timescales as a function of the chirp mass of 
merging systems. The dashed horizontal and vertical lines 
indicate threshold values, Tm.thrc and A^thro, so that 70% of 
the systems have rm > Tm.thro and 75% have A4 < Mthre- 



merging pairs at the time of formation of the compact bi- 
nary systems. The inspiral part of the adopted waveform 
depends on the BH masses and it is strictly applicable 
only to (BH-BH) binaries in quasi-circular orbit. From 
Fig.[4]wc sec that only 10% of the systems have e > 0.15. 
Furthermore, using the relation between the orbital sep- 
aration and the eccentricity due to gravitational wave 
emission [13, Ell, we find that by the time (BH-BH) bina- 
ries enter the Advanced LIGO and ET bands (10 Hz and 
1 Hz, respectively), the eccentricities are smaller than 
10~^ for all systems, and consequently the assumption of 
quasi-circular orbits is well justified. 



IV. BLACK HOLE BINARIES AS GW SOURCES 

To evaluate the background produced by a cosmolog- 
ical population of coalescing (BH-BH) binaries, we use 
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FIG. 4: Eccentricity as a function of the black hole mass ratio 
for merging systems. 



the family of waveforms obtained in [l5|, which model 
the inspiral, merger and ring-down phases for the coa- 
lescence of non spinning black holes. These waveforms 
refer to the leading harmonic of the gravitational signal 
{i = 2,m = ±2), which is the dominant contribution for 
low mass ratios (qbh !)• In the frequency domain, the 
signal has the form 



M/)=^cff(/)e'*-(^\ 



k 


0>k 


h, 







6.6389 X 10^' 


-1.0321 X 10"' 


1.0979 X 10"' 


1 


1.3278 


-2.0642 X 10"' 


2.1957 X 10"' 


2 


1.1383 


-1.7700 X 10"' 


4.6834 X 10"^ 


3 


1.7086 


-2.6592 X 10"' 


2.8236 X 10"' 



TABLE II: Values of the constants which appear in the wave 
amplitude (eq. [2]) taken from [isj . 
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FIG. 5: (color online) Single source signal as a function of 
frequency for a (BH-BH) binary with total mass M — 20Mq, 
symmetric mass ratio -q = 0.25, located at a distance of 
10 Mpc. 



where / is the emission frequency. The wave amplitude 
is given by 



wC 



-7/6 
-2/3 



if / ^ Iraerg 

if fmerg ^ f ^ fring 
if fring ^ I ^ I cut 



(2) 

where the constants (fmerg, frmg: f cut): which identify 
the frequency regions where the emitting system is in- 
spiralling, merging and oscillating, are 



Jmerg — _ , , ('J ) 

fring 
fcut 





+ boV^ 


\- Co 




ttM 




airf 


+ foir/ H 






7:M 










ttM 



In these expressions 77 = Mi^bhA'^2,bh/-M^ is the sym- 
metric mass ratio, and the coefficients a^, 6^, and (with 
/c = 0, 1, 2, 3) are given in Table HIl The constants C and 
If, and the function C, which characterize the wave am- 



plitude, are 



C 



w 



^_2/3 fV6 V24 

u n Jmerg 



fr 



ing 



7Ta 

^ \ fmerg , 
1 \ cr 



24 

-2/3 



1/2 



(4) 



27r; (/-/™<,)2+a2/4' 



where cr = (02??^ -I- 62'? + C2)/ttM. We do not need to 
compute the phase of the signal since the single source 
GW spectrum depends on the squared GW amplitude 
(see eq. Ilip . The model assumes optimal orientation of 
the detector with respect to the emitting source. 

In Fig. [5] we plot the dimensionless amplitude fh{f) 
as a function of frequency, for a system with total mass 
M = 2OM0, symmetric mass ratio r] ~ 0.25, located at a 
distance of 10 Mpc. The frequency limits, which identify 
the three regimes (inspiral, merger and ring-down), are: 

810 Hz, and fcut = 1042 Hz 



f. 



merg 



405 Hz, /, 



ring 



(points on the curve). These frequencies are inversely 
proportional to the total mass of the binary. Hence, for 
the most massive binaries in the simulated sample, with 
M = 4OM0 (see Fig. d]), f^nerg, frzng, and fcut are a 
factor of 2 smaller than the values plotted in the figure. 
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V. FROM STAR FORMATION TO BINARY 
FORMATION RATE 

In this section we derive the evolution of the birth 
rate of binary systems from the comoving star forma- 
tion rate density as a function of redshift, inferred from 
the simulations in [l9l |. These cosmological simulations 
are characterized by an improved treatment of metal en- 
richment and the stellar IMF is assigned depending on 
the gas mctallicity. In particular, Population II/I stars 
form in the mass range [0.1 — IOOJA/q according to a 
Salpetcr IMF in regions which have been already pol- 
luted by the first metals and dust grains to a metallicity 
Z> Z„ = [10-*^ - W)-^\Zq [H-ieii- Below this thresh- 
old, gas cooling is inefficient and the star formation pro- 
cess favors the formation of very massive (Population III) 
stars, characterized by a top-heavy IMF (for more details 
on the numerical scheme we refer to [l9|). In this work, 
we are only interested to Population II/I stellar progen- 
itors of black holes. 

Starting from the star formation rate density at a given 
2, Pi,(z) (expressed in units of Mq yr~^Mpc~'^), we derive 
the binary birth rate per comoving volume (expressed in 
units of yr"^ Mpc~^) as. 



i?bin(z) 



dR 
dtdV 



{z)^^x^x U^, (5) 
(to*) 2 



where /bin is the binarity fraction (which we take to be 1), 
(to*) is the average stellar mass and fsim is the fraction 
of binaries simulated by the population synthesis code 
ScBa. The latter quantity accounts for the fact that, 
while in the original simulation of stars are assumed 
to have masses in the range [O.l-IOOJM©, in SeBa we 
initialize only binary systems with primary mass in the 
range [8-100]Mq, in order to increase the statistics on 
double (BH-BH) binaries. Thus, the fraction of simulated 
systems is, 



_ /8°°rfA/prim$(Mp,i„,) 
/sim -- 100 r ^/ 7. r \ ' 



(6) 



and the average stellar mass is, 



i$(Mpnm) 



/o°°dAfp,i,n*(M, 



(7) 



prim / 



Following [H , we assume that a Z AMS binary forms at a 
redshift z^; after a time interval the system has evolved 
into a (BH-BH) binary. Consequently, the redshift of 
formation of the degenerate binary system, zj is defined 
through t{zf) = t{zs) + Tg. Once the (BH-BH) binary 
system is formed, it evolves according to gravitational 
wave emission until, after a time interval t^, it eventually 
coalesces. The redshift Zc at which coalescence occurs is 
defined by t(zc) = t(zf) + r^. The number of (BH-BH) 
systems formed per unit time and comoving volume at 




FIG. 6: (color online) Redshift evolution of the total binary 
birth rate (solid red line) and of (BH-BH) birth and merger 
rates (dotted and dashed, respectively). 



redshift Zf is, 



pbirth \ 
^(BH-BH) \Zf) 



dT„ 



tizf)-t{zp) 



dr. 



(8) 



^(BH-BH) i?bin(^s) , x 

X r-]3(BH-BH)(T-s,Tmj, 



N 



{1 + z,)' 



where zp defines the onset of star formation (zp ^ 20 
in the simulation), A'^(bh-bh) is the number of (BH-BH) 
systems, N is the total number of simulated binaries, and 
P {BH-BH) {TsiTm) is the joint probability distribution of 
delay times. Similarly, the number of (BH-BH) systems 
per unit time and comoving volume which merge at red- 
shift Zc is. 



Amorgor / x 
^'■(BH-BH)V'^cy 



t{Za)-t{zF) 



dT„ 



t(z^)-T^-t(zF) 



^(BH-BH) i?bin(^s) , . 

X — tP(BH~BH) [Ts,Trn)- 



N 



(l + Zs)' 



dTs 

(9) 



In Fig. ini we show the redshift evolution of the binary 
birth rate (solid line). Of all these systems, only 1.7 % 
form a (BH-BH) binary. The evolution of the (BH-BH) 
birth and merger rates is also shown in the figure (dashed 
and dotted lines, respectively). Since Tg is relatively short 
(less than ~ 6 Myr) , the evolution of -RJbh-bh) simply 
a scaled-down version of -Rbin- Conversely, there is a 
shift in the evolution of -R^b'^h^bh) ^lii^h is due to the 
long merger timescales, Tm > [1 — 5] Gyr. 

For the sake of comparison, we calculate the Galactic 
birth and merger rates for (BH-BH), (BH-NS) and (NS- 
NS) systems extracted from the simulation of Model A. 
To compute these rates, we normalize the total number 
of core-collapse SNe that we find in the simulation to an 
estimated Galactic supernova rate of 1 x 10~^yr~^ [6^ . 
The resulting values are presented in Table Hill and are in 
good agreement with the current literature on the subject 
(see [3. M [26l. l63j and [S^ and references therein). 
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Galactic Birth/Merger Rates 


Type 


Birth rates (yr ^) 


Merger Rates (yr ^) 


(NS-NS) 


8.2x10"^ 


2.0x10"'' 


(BH-NS) 


5.3x10"^ 


6.2x10"" 


(BH-BH) 


9.5x10"^ 


1.8x10"" 



TABLE III: Galactic Birth/Merger rates obtained from the 
simulation of Model A normalizing to a Galactic supernova 
rate of 1 x 10"^yr~^ (see text). 



VI. GWB FROM BLACK HOLE BINARIES 

Following 0, [l3] , we write the spectral energy density 
of the GWB produced by a population of (BH-BH) bi- 
naries as, 



dE 



dSdfdt 



c^^(bh'^bh)( 



dE 
dSdf 



(10) 



where ^^^(bh^bh) = ^(bh-bh) and the locally mea- 
sured average GW energy flux from a single source at 
rcdshift z is 



^dSdf' G2-' ^ 



zf\h[f{l + z) 



(11) 



Here / = fe{l + z)^^ is the redshifted emission fre- 
quency and h is the amplitude of the GW signal given in 
Eq. ([T]). The GWB is conventionally characterized by the 
dimensionless quantity ^Icwif) ^ Pcr~^{dpgw / d\og f), 
which is related to the spectral energy density by the 
equation, 



f 



C?Pc 



dE 



dSdfdt 



(12) 



where p^^ = SH^/SnG is the cosmic critical density. 

In Fig. [7] we plot I^gw, as a function of the obser- 
vational frequency, for the reference model A. The cu- 
mulative signal is the result of the emission during the 
inspiral (dashed line), merger (solid line) and ring-down 
(dotted line) phases of the coalescence processes. In the 
frequency range 10 Hz < / < 200 Hz. the signal is 
dominated by the inspiral phase which reaches a max- 
imum amplitude of Oqw = 7.8 x 10"^° at a frequency 
of ^ 200 Hz. Above this limit, a further increase in the 
signal is driven by the emission during the merger phase, 
which reaches a maximum J^gw = 2.1 x 10~^ at 540 Hz. 
This is not surprising, since a significant portion of GW 
energy is radiated during the merger. At larger frequen- 
cies, the signal drops with a minor contribution coming 
from the ring-down phase which follows the final coales- 
cence of the two black holes in each binary. 

In Fig. [3 the GWB signal is compared with the 
foreseen sensitivity curves for advanced LIGO/ Virgo 
(ALIGO) and for the Einstein Telescope with two dif- 
ferent design configurations (ET-B and ET-C) and as- 
suming 1 yr or 3 yrs of integration time. In particular. 
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FIG. 7: (color online) The closure energy density, Qqw, gen- 
erated by (BH-BH) coalescing binaries in model A, plotted 
as a function of the observational frequency. The three con- 
tributions coming from the inspiral (dashed), merger (solid) 
and ring-down (dotted) phases are plotted separately. The 
black solid line is the total GWB signal. The shaded regions 
in the three panels indicate the foreseen sensitivities for dif- 
ferent interferometers and integration times (see labels on the 
figures) . 



ALIGO with 3 years of integration might sample only 
a small portion of the inspiral phase with (S/N) ratios 
below the detection threshold (see section I VI A[) . The 
full inspiral and merger phases might be observed with 
ET, even with an integration time of 1 year. The best 
configuration appears to be ET-B (see discussion below) 
which amplifies the sensitivity at larger frequencies. In 
section IVl Ai we estimate the detectability of the signal 
in a quantitative way and in section IVIBI we discuss the 
dependence of the GWB and of its detectability on some 
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key physical parameters. 



A. Detectability 

We consider the design sensitivities of second gener- 
ation interferometric detectors, Advanced LIGO/ Virgo 
(rij , in a configuration of zero-detuning of the signal re- 
cycling mirror, with high laser power. For the third gen- 
eration interferometer ET, we consider two target sensi- 
tivities. The first configuration, ET-B, is an underground 
based design, incorporating long suspension, cryogenics, 
and signal power recycling. The second configuration, 
ET-C, is called Xylophone configuration and merges the 
output of two detectors specialized in different frequency 
bands (for more details see [63. Issj ) . 

It is known that the detection strategy for continuous 
GWB signals is to cross-correlate the output of two de- 
tectors that are assumed to be sufficiently well separated 
that their noise sources are largely uncorrelated [T^l ■ 

The statistical nature of the background depends on 
the duration of the signal, on the event rate and on the 
lower frequency bound of the detector /l , the so-called 
"seismic wall" (see for a discussion [6a])- It has recently 
been suggested that, whenever the signal is very short 
(as it is expected for binaries at high frequencies), and of 
sufficiently high amplitude, the resolution of the detector 
or of the data analysis method will enable to individually 
identify the signals and subtract it from the data (67j . 

Assuming that for ALIGO /l = 10 Hz, the background 
predicted for the reference model A would be character- 
ized by a duty cycle of 5 x 10~^. For ET, which has 
a lower /l = 1 Hz, the duty cycle would be ^ 2. Ac- 
cording to [il] , the cross-correlation method is found to 
be nearly optimal for duty cycles > 10"'^. Therefore, in 
what follows we estimate the detectability of the signal 
using this method, as in Q , and we do not distinguish be- 
tween the resolvable and the non-resolvable components 
1671. We use the cross correlation statistics to calculate 



9^ 
507r4 



df 



(13) 



(S/N) 


Case-II [Case-I] 


GWB 


ALIGO 


ET-B 


ET-B 


ET-C 




(3 yr) 


(3 yr) 


(lyr) 


(lyr) 


coalescence 


1.9[0.7] 


316[118] 


182 [68] 


275[103] 


inspiral 


1.9[0.7] 


310[116] 


179[67] 


274[103] 


merger 


7.5xl0"''[6.4xl0"*] 


26 [9] 


15[5.6] 


5[1.9] 


ring-down 


8.1xl0"-'[3.3xl0"''] 


1.5[0.6] 


0.9[0.3] 


0.4[0.1] 



TABLE IV: The (S/N) ratio for second and third generation 
detectors, assuming different integration times and detector 
separation/orientation (see text). The values refer to model 
A (see table [I]) and have been computed considering the cu- 
mulative signal (coalescence) and separate contributions from 
the inspiral, merger and ring-down phases. 

the optimized S/N for an integration time T as given by 



where Pi{f) and ^2(/) are the power spectral noise den- 
sities of the two detectors and 7 is the normalized overlap 
reduction function, which quantifies the loss of sensitiv- 
ity due to the separation and the relative orientation of 
the detectors. 

We have computed the (S/N) ratio assuming differ- 
ent integration times (1-3 years) and detector separa- 
tion/orientation (casc-I and -II). For ALIGO, case-I con- 
siders the LIGO Hanford/Livingston pair using 7 in the 
form given by Eq. (3.26) in [gij. For ET-B and ET-C, 
case-I adopts the constant value 7 = —3/8 which applies 
to two ET detectors operating in the frequency range 
[1-1000] Hz [6J]. Case-II is the same for all detectors 
and represents a pair of aligned equivalent detectors sit- 
uated within several km. This optimal case corresponds 
to 7 = 1. 

The resulting (S/N) ratios are reported in Table IIVI 
The values refer to model A (see table |T]) , and are ob- 
tained assuming a threshold signal-to-noise ratio of 3 
which corresponds to a false alarm rate of 10% and to 
a detection rate of 90% (for more details see eq. 19 in 

El)- 

The highest {S/N) ratios are obtained with optimal 
orientation (case-II) and longer integration times. These 
conditions would not allow ALIGO to detect the coales- 
cence signal. The increase in sensitivity foreseen for ET 
would enable the detectability of two portions of the sig- 
nal, the inspiral and the merger phases, with [S/N] > 5, 
independently of the adopted separation/orientation and 
integration times. 



B. Dependence on physical parameters 

In this section we analyze the dependence of the GWB 
on some key physical parameters that affect the (BH- 
BH) birth/merger rates, i.e. the adopted kick velocity 
distribution, the CE parameter and core mass threshold 
for black hole formation, mthre.BH- 

In Table IVl we list a set of models (first column) which 
differ from the reference model A by the variation of a 
single parameter, indicated in the second column. For 
each model, in column 3 and 4 we tabulate the Galactic 
(BH-BH) birth and merger rates. Similarly to what done 
in Section |V] (Table IIIip , these values are obtained by 
normalizing the supernova rate predicted by each SeBa 
run to the Galactic rate of 1 x 10~^ yr~^. In column 5 
we give the local merger rate, which is the z — value of 
the cosmic merger rate (see for instance Fig. [6] for model 
A) obtained as described in Section |Vl Eq. (0). 

In Table 8 of ref . [t^] , the Galactic rates obtained with 
different population synthesis codes, based on different 
assumptions on stellar/binary evolution, are compared. 
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EXPLORATION OF THE PARAMETER SPACE 


Model 


Modified Parameter 


Galactic (BH-BH) 
BR (Myr^^) 


Galactic (BH-BH) 
MR (Myr"^) 


Local (BH-BH) 
MR (Mpc~^Myr"') 


A 


reference model 


95 


1.8 


5.3 X 10"^ 


Bl 


acE^=0.5 


95 


16 


2.9 X 10"' 


B5 


acE^=4: 


96 


1.2 X 10"' 


3.3 X 10"^ 


CI 


Paczynsky distribution (cr — 150 km s~') 


120 


1.7 


5.1 X 10"^ 


C2 


Paczyfisky distribution (cr — 600 km s"') 


71 


2.4 


6.6 X 10"^ 


C5 


Maxwellian distribution (cr — 200 km s ) 


88 


1.8 


5.3 X 10"^ 


C7 


Maxwellian distribution (cr — 600 km s"') 


53 


3.2 


8.8 X 10"^ 


C8 


No kick distribution 


220 


1.7 


5.3 X lO-"" 


Fl 


mthrc.BH = 8.5M0 


140 


7.6 


2.0 X 10"' 


F2 


mthrc.BH = 7.6M0 


160 


15 


3.6 X 10"' 


F3 


mthrcBH = 5.5M0 


240 


40 


8.5 X 10"' 



TABLE V: In columns 3 and 4 we give the Galactic Birth (BR)/Merger rates (MR) for the different models, obtained by 
varying one key parameter (indicated in colum 2) with respect to the reference model A (see text). In column 5 we give the 
local merger rate, which is the z = value of the merger rate (see for instance Fig. [6] for model A) obtained as described in 
Section [Vl Eq. ^. 



The values reported in that table range between 0.01 and 
250 Myr~^, and values that are considered as "realistic" 
range within 0.01 and 20 Myr~^. If we compare these 
estimates with the Galactic merger rates given in column 
4 of Table |Vl we see that our values are in the realistic 
range, except for the model F3, which we will discuss 
below. 

We first consider models where we vary the common 
envelope parameter (models Bl and B5). This parame- 
ter does not affect the birth rate but controls the num- 
ber of merging systems: in fact, a larger (smaller) CE 
parameter, such as in model B5 (Bl), generates binaries 
which, at the end of the CE phase, are characterized by 
larger (smaller) orbital separations and therefore longer 
(shorter) merger timescales. The resulting fraction of 
merging systems in model B5 is only 0.12% of the total 
(BH-BH) binaries, more than a factor 10 smaller then for 
the reference model A. 

The adopted shape and velocity dispersion of the kick 
distribution affect both the birth and merger rates (see 
models CI, C2, C5 and C7). In the first two models, we 
assume a Paczynsky distribution for kick velocities, 

P{u)du = - ^" (14) 

TT (1 + W^)^ 

where u = v/a and the dispersion velocity a is, respec- 
tively, 150 and 600 km s~^ for models CI and C2. In 
models C5 and C7 we consider a Maxwellian distribu- 
tion, 

P{v)dv = J-v^e'^ (15) 

V TT 

with cr = 200 and 600 km s"^, respectively. For both 
distributions, the birth rate decreases with increasing a. 
However, the fraction of formed binaries which coalesces 



increases slightly with a. The first effect is a direct con- 
sequence of the disruption of the binary after the SN ex- 
plosion (it is more likely that a strong kick disrupts the 
system). The increase of the merger rate with a is due to 
the net effect of kicks on the orbit of systems which are 
not disrupted: while the semi-major axis is left mostly 
unchanged, the eccentricity is greatly enhanced, favour- 
ing the orbital decay due to GW emission. As a result, 
in model C2 the number of merging pairs is a factor 2 
larger than in model A (3.4% of the total sample). In the 
extreme case of null kicks, such as in model C8, we find 
a birth rate which is a factor 2.3 larger than in model A. 
Yet, the corresponding merger rates are comparable in 
the two models. 

Finally, we discuss the dependence on the core mass 
threshold for BH formation (models F1-F3). This pa- 
rameter appears to be the most important one for the 
GWB. In fact, a reduction in wthrcBH from the reference 
value of IOMq (model A) to 5.5Mq (model F3) leads to 
a birth rate which is a factor 2.5 larger. The amplifica- 
tion in the merger rate is even more dramatic, by more 
than a factor 20. This is due to the evolutionary path 
followed by massive stellar progenitors. In particular, a 
smaller mthrc.BH allows the formation of BHs from lighter 
progenitor stars; the latter experience a smaller amount 
of mass loss (see Section |IT|) forming close binary pairs 
characterized by shorter merger timescales. 

The largest differences with respect to model A are 
found for models B5 and F3 which provide a sort of lower 
and upper limits to the (BH-BH) merger rate and GWB. 
In Fig.|n]we show the predicted redshift evolution of (BH- 
BH) merger rate in the two models as compared to model 
A. The upper limit labelled as 'LIGO S5' shows the con- 
straint derived in [l[ from approximately 2 years of LIGO 
data (run S5) on the merger rate of systems with compo- 
nent masses in the range IQMq — 28Mo. Even the most 
optimistic model F3 predicts a local merger rate which is 
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FIG. 8: (color onlme)f2Gw generated during the inspiral (up- 
per panel), merger (centra/ panel) and ring-down (fo«;er panel) 
phases in models F3 (dashed), B5 (dotted) and A (solid). In 
all panels, the two shaded regions indicate the foreseen sensi- 
tivities of ALIGO and ET-B assuming 3 years of integration. 



more than a factor 2 smaller than the upper limit inferred 
by the data analysis of the LIGO S5 run. 

The longer (shorter) merger timescales predicted in 
model B5 (F3) lead to a reduction (amplification) of the 
overall cosmic merger rate, shifting it to smaller (larger) 
redshifts. This, in turn, affects the amplitude and fre- 
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FIG. 9: (color online) Redshift evolution of (BH-BH) merger 
rates for models F3 (solid line), A (dotted line) and B5 
(dashed line). The upper limit labeled 'LIGO S5' shows the 
constraint recently derived in [J] (see text). 



quency range of the resulting GWB spectra, as shown in 
Figs, il and [101 

In Fig. ISl wc plot the contributions to ficw generated 
during the inspiral (upper panel), merger (central panel) 
and ring-dovirn (lovifer panel) phases comparing models 
F3, B5 and A. As expected, model F3 generates the 
strongest signals. It is also evident that the merger and 
ring-down signals in model F3 extend to lower frequen- 
cies with respect to models A and B5. This is due to 
the shorter merger timescales which allow a larger num- 
ber of (BH-BH) binaries to reach the final coalescence at 
larger redshifts (see Fig. [9|) , emitting signals contribut- 
ing at smaller observational frequencies. Similarly, the 
differences among models A and B5 can be traced back 
to the merger timescales which, for models B5, confines 
(BH-BH) coalescence to redshifts z < 1. It is also in- 
teresting to note that while models A and B5 show a 
similar behaviour at the largest frequencies, model F3 
systematically extends to larger frequencies. In fact, the 
smaller black hole masses predicted in model F3 (as a 
consequence of the smaller TOthrc.Bn) lead to larger /merge, 

/ring! and /eut- 

The shaded region in Fig. [TU] illustrates the largest vari- 
ations of the GWB among the models, and can be viewed 
as an indication of the uncertainty affecting its estimate. 
For the models we consider, the peak amplitude in the 
closure energy density ranges within 10~^° < f^cw < 
5 X 10"*^ at frequencies 470 Hz < / < 510 Hz. 

Table IVII quantifies these differences in terms of the 
predicted signal-to-noise ratio. For the sake of compari- 
son, we consider the same integration times and detector 
configurations/orientations as in Table HVl which refer to 
model A. The major difference with respect to the data 
in Table HVl is that, for model F3, ALIGO has a chance to 
probe the inspiral part of the GWB with an integration 
time of 1 yr and (S/N) larger than 10. 
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(S/N) 


Case-II [Case-I] 


O W-D 


iVLoaei 






til 1 - Jj (^i yr ) 


rj L-D [6 yr ) 


inspiral 


F3 


29[10] 


50[19] 


4720[1770] 


8176[3065] 


inspiral 


B5 


6.4xl0~^[2.4xl0"^] 


0.1[4.1xl0"^] 


10[3.9] 


18[6.7] 


merger 


F3 


l.l[1.3xl0"^] 


1.9[2.3xl0"^] 


442[165] 


766[287] 


merger 


B5 


S.OxlO--' [3.3x10"''] 


5.2x10-' [5.6x10"''] 


1.2[0.4] 


2.1[0.8] 


ring-down 


F3 


O.liT.OxlO-*] 


0.2[1.2xl0"''] 


32[12] 


56[21] 


ring-down 


B5 


3.3xl0"*[1.7xl0"''] 


5.7xl0"*[2.9xl0"''] 


8.2xl0"^[3.0xl0"^] 


0.1[5.3xl0"^] 



TABLE VI: The (S/N) ratio for second and third generation detectors assuming different integration times and detector 
separation/orientation. These values have been computed considering separate contributions from the inspiral, merger and 
ring-down phases for models F3 and B5. The numbers in square brackets refer to a non optimal orientation of the detectors 
(see the discussion in Sec. lVI A)l . 
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FIG. 10: (color online) The closure energy density, ficw, for 
(BH-BH) binaries predicted by models F3 (dashed line) and 
B5 (dotted line) is compared to that of model A (solid line). 
The shaded region can be considered a measure of the un- 
certainty on ficw. The two black solid lines indicate the 
foreseen sensitivities of ALIGO and ET-B assuming 3 years 
of integration. 



VII. DISCUSSION AND CONCLUSIONS 

Many astrophysical processes which control the for- 
mation and evolution of (BH-BH) binaries, starting from 
their stellar progenitors, are still poorly understood. Of 
particular relevance in this respect, are the amount of 
mass loss by massive stars, pre-supernova and supernova 
evolution, the effects of mass transfer among the two 
companion stars on the subsequent evolution of the sys- 
tem. In the present study, we have considered a reference 
model (A), which adopts a set of "standard" conserva- 
tive assumptions, reproducing the observed prop erties of 
single Wolf-Rayet stars and double pulsars ( [3^ and ref- 
erences therein). Furthermore, we have explored a wide 
range of parameters on which our simulations depend 
(see Table 1), to extract those that have the largest im- 
pact on the GWB, i.e. the common envelope parameter, 
the core mass threshold for BH formation and the kick 



velocity distribution. Varying these parameters, we iden- 
tify two models, B5 and F3, which produce, respectively, 
the smallest and largest gravitational wave background. 

When normalized to a Galactic star formation rate, 
the Galactic birth and merger rates computed for the 
considered models are in goo d ag reement with the results 
of independent studies (see [70| and references therein, 
and the discussion in Sec. |V|) . 

To model the sin gle source emission, we use the hybrid 
waveforms given in [15| , which refer to non-spinning (BH- 
BH) binaries. This model has been improved in [7l| with 
the inclusion of the effect of non-precessing spins, and of 
a more accurate modeling of the non-spinning case. A 
comparison of the GWB obtained with the waveforms of 
[Tit and of [7l|, for our reference model A, assuming the 
single spin parameter x — (non spinning case), shows 
no significant differences. Even assuming that all (BH- 
BH) binaries have x = 0.85 (which is the extreme case 
considered in [7l| ) , the GWB shows some difference only 
above ~ 1 kHz, in a region where even ET could not de- 
tect it. This evidence has been reported also in In 



71] it is stated that the non-precessing waveforms are ef- 
fectual in capturing also precessing binaries, expecially in 
the comparable-mass regime. Since most of our systems 
have g < 1.1 (see figure^]) the inclusion of spin precession 
is not expected to significantly change our results. 

As mentioned in the introduction, the same waveforms 
have been recently used by Q to estimate the gravi- 
tational wave background. This study adopts a fixed 
merger time distribution function, an average chirp mass 
to describe all the (BH-BH) systems, and normalize the 
binary merger rate to a "local" merger rate (LMR); 
the latter quantity is computed multiplying the Galac- 
tic merg er rates predicted by Population Synthesis Mod- 
els, 



51 



by the average number density of Milky Way- 
type galaxies in the Local Universe (assumed to be 0.01 
Mpc""^). Thus, there are many differences between this 
study and the present one. As explained in section |Vl 
we do not normalize the cosmic merger rate to a "re- 
scaled" Galactic rate; indeed, we compute it integrat- 
ing the star formation history over the birth and merger 
time distribution functions predicted by each SeBa run. 
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The resulting local rates, at 2; = 0, are shown in col- 
umn 5 of Table |Vl and are a factor 2-3 (depending on 
the model) larger than what would be obtained from the 
corresponding Galactic merger rate (column 4), applying 
the procedure of In fact, due to their long merger 
times, a large fraction of the (BH-BH) progenitors form 
at z > 1 close to the peak of the cosmic star formation 
rate. In addition, our analysis shows that there exist 
correlations between the distributions of merger times, 
hence the merger rate, and the distribution of the chirp 
masses (see Fig. [3]). Thus, the points of the detectable 
parameter space explored by @ , identified by an average 
chirp mass and a LMR, do not have the same probabil- 
ity to represent a physical model. Our study shows that 
variations of key physical parameters produce correlated 
effects on the merger time and chirp mass distributions 
and; as a consequence, on the amplitude and the spectral 
energy distribution of the gravitational wave background 
(see, for instance, the discussion on model F3 in section 
IVIBp . These properties can be appreciated only using 
the full rich information provided by Population Synthe- 
sis Models. 

Our main results can be summarized as follows: 

• For the reference model A, the sample of simulated 
(BH-BH) binaries is characterized by BH masses 
which vary between ~ 6Mq and 20Mq, with 
the largest concentration in the range [10 — ISJA/q. 
The formation of (BH-BH) binaries from their stel- 
lar progenitors is characterized by relatively short 
timescales, ~ 3.5 — 6 Myr and by a wide interval of 
semi-major axes ranging between lOi?0 to sev- 
eral thousands of Rq. Only 2% of the formed (BH- 
BH) binaries are able to merge within the Hubble 
time. These systems are characterized by semi- 
major axes < 2Oi?0, and black hole mass ratios 
close to 1. The majority of these systems (70%) 
have merger timescales > 6 Gyr. As a result, (BH- 
BH) Galactic birth and merger rates are, respec- 
tively, 9.5 X 10"^ yr-i and 1.8 x 10"^ yr'^. On 
cosmic scales, the (BH-BH) birth rate closely fol- 
lows the shape of the cosmic star formation rate 
(although with a significantly reduced amplitude); 
conversely, the merger rate shows a significant time 
delay, and it is negligible beyond z ^ 2. The above 
conclusions mostly depend on the adopted common 
envelope parameter and core mass threshold for BH 
formation [acE^ = 2 and mthrc.BH = lOil/© in 
model A). 

• An increase of the CE parameter to acE^ — 4, 
as in model B5, generates (BH-BH) binaries with 
larger orbital separation, reducing the (BH-BH) 
merger rate (by a factor of 10 for the Galactic 
value) and confining the mergers to occur at z < 1. 
Conversely, a reduction in mthrc.BH to 5.5M0, as 
in model F3, leads to an increase of the Galactic 
merger rate by more than a factor 20. Variations 
in these physical parameters also affect the distri- 



bution of black hole masses. 

• The GWB is characterized by a peak amplitude in 
the range 10" < flow < 5 x 10~^ at frequencies 
470 Hz < / < 510 Hz, when the uncertainties on 
some key physical parameters are considered (see 
Figlini). 

• Advanced LIGO /Virgo have a chance to detect the 
GWB from the inspiral only in model F3, which 
predicts the highest merger rate; third generation 
detectors like ET, would detect the inspiral GWB 
with high (S/N) for models spanning the region 
from model A to model F3 in Fig. [10] (see Tables ITVl 
andim. 

• The merger contribution to the GWB could be de- 
tected only by ET. From Fig. [Hand Tables HVl and 
IVIi we see that ET-B could detect this contribu- 
tion for models spanning the region from model A 
iiS/N) > 6) to model F3 {{S/N) > 165), with 1 
yr integration. 

• The ring-down contribution could be detected by 
ET-B with 1 yr integration only for model F3 
((^/iV)> 12). 

We find that the amplitude of the GWB is very sensi- 
tive to the adopted core mass threshold for BH formation 
(models Fl to F3). This opens up the possibility to con- 
strain the uncertain physics related to the final stages of 
the evolution of massive stars using observational con- 
straints on the associated gravitational wave emission. 

Finally, we would like to mention that, according to 
a recent study [5lj, (BH-BH) binaries formed in a low- 
metallicity environment are characterized by higher coa- 
lescence rates and chirp masses than their solar metallic- 
ity counterparts. Indeed, at sub-solar metallicities stars 
arc more compact (smaller radii), experience reduced 
mass loss (larger masses) and BHs can form by direct 
collapse of the progenitor (no kick due to the SN explo- 
sion) . These changes result in an increase of the galactic 
and local merger rate. Metallicity-dependent evolution- 
ary tracks are currently being implemented in SeBa [t^ , 
and their effects on the gravitational wave background 
will be considered in a future study. 
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